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Abstract 

We extend our formulation of relativistic three-nucleon Faddeev equations to include both pairwise in- 
teractions and a three-nucleon force. Exact Poincare invariance is realized by adding interactions to the 
mass Casimir operator (rest Hamiltonian) of the non-interacting system without changing the spin Casimir 
operator. This is achieved by using interactions defined by rotationally invariant kernels that are functions 
of internal momentum variables and single-particle spins that undergo identical Wigner rotations. To solve 
the resulting equations one needs matrix elements of the three-nucleon force with these properties in a 
momentum-space partial-wave basis. We present two methods to calculate matrix elements of three-nucleon 
forces with these properties. For a number of examples we show that at higher energies, where effects of rel- 
ativity and of three-nucleon forces are non- negligible, a consistent treatment of both is required to properly 
analyze the data. 

PACS numbers: 21.45.-v, 21.45.Ff, 25.10.-Hs, 24.10. Jv 
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I. INTRODUCTION 

High precision nucleon-nucleon potentials such as AV18 CDBonn Nijm I, II and 93 
provide a very good description of the nucleon-nucleon data set up to about 350 MeV. When 
these forces are used to predict binding energies of three-nucleon systems they underestimate the 
experimental bindings of and ^He by about 0.5-1 MeV [4, 5|. This missing binding energy can 
be restored by introducing a three-nucleon force into the nuclear Hamiltonian ^. 

Also the study of elastic nucleon-deuteron scattering and nucleon induced deuteron breakup 
revealed a number of cases where the nonrelativistic description using only pairwise forces is in- 
sufficient to explain the data. Generally, the studied discrepancies between a theory using only 
nucleon-nucleon potentials and experiment become larger with increasing energy of the three- 
nucleon system. Adding a three-nucleon force to the pairwise interactions leads in some cases to a 



better description of the data. The elastic nucleon-deuteron angu 



ar distribution in the region of 



3,Ei. 



The clear discrepancy in 



its minimum and at backward angles is the best studied example 
these angular regions at energies up to w 100 MeV nucleon lab energy between a theory using only 
nucleon-nucleon potentials and the cross section data can be removed by adding a modern three- 
nucleon force to the nuclear Hamiltonian. Such a three-nucleon force must be adjusted with each 
nucleon-nucleon potential separately to the experimental binding of "^H and ^He 6-^. At energies 
higher than ~ 100 MeV current three-nucleon forces only partially improve the description of cross 
section data and the remaining discrepancies, which increase with energy, indicate the possibility 
of relativistic effects. The need for a relativistic description of three-nucleon scattering was also 
raised when precise measurements of the total cross section for neutron-deuteron scattering ^ were 
analyzed within the framework of nonrelativistic Faddeev calculations [10]. Nucleon-nucleon forces 
alone were insufficient to describe the data above ~ 100 MeV. The effects due to relativistic kine- 
matics considered in were comparable at higher energies to the effects due to three-nucleon 
forces. These results showed the importance of a study taking relativistic effects in the three 
nucleon continuum into account. 
In 
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12l | the ffist results on relativistic effects in the three-nucleon continuum have been 
presented. The dynamics was defined by a three-nucleon center of momentum Hamiltonian or 
mass operator including only pairwise interactions. The mass operator was used to calculate three- 
nucleon scattering observables. The input to that approach is a "Lorentz boosted" nucleon-nucleon 
potential, which generates the nucleon-nucleon t-matrix in a moving frame by solving a standard 
Lippmann-Schwinger equation. To get the nucleon-nucleon potential in an arbitrary moving frame 



2 



one needs the interaction in the two-nucleon center of momentum system, which appears in the rela- 
tivistic nucleon-nucleon Schrodinger or Lippmann-Schwinger equation. The relativistic Schrodinger 
equation in the two-nucleon center of momentum system differs from the nonrelativistic Schrodinger 
equation just by the relativistic form for the kinetic energy. Current realistic nucleon-nucleon po- 
tentials are defined and fit by comparing the solution of the nonrelativistic Schrodinger equation 
to experimental data. Up to now nucleon-nucleon potentials refitted with the same accuracy in 
the framework of the relativistic nucleon-nucleon Schrodinger equation do not exist. Such refitting 
can be, however, avoided by solving a quadratic integral equation whose solution is a relativistic 
potential which is phase-equivalent to a given input high-precision nonrelativistic nucleon-nucleon 
potential [l^. An alternative equivalent approach towards a relativistic nucleon-nucleon t-matrix 



in another frame is provided in 



In our previous studies with only nucleon-nucleon interactions we found that when the non- 
relativistic form of the kinetic energy is replaced by the relativistic one and a proper treatment of 
the relativistic dynamics is included, the elastic scattering cross section is only slightly influenced 
by relativity. Only at backward angles and higher energies are the elastic cross sections increased 
by relativity [llj . It is exactly the region of angles and energies where the effects of three- nucleon 
forces are also significant [8]. Also, for nucleon-deuteron breakup reactions regions of phase space 
were found at higher energies of the incoming neutron where relativity significantly changes the 



breakup cross sections 
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Ig ]. For some spin observables large effects due to relativity and three- 



nucleon forces have been reported in nucleon-deuteron breakup for an incoming deuteron energy of 
270 MeV, some of which seem to be supported by proton-deuteron data [17]. These observations call 
for three- nucleon continuum relativistic Faddeev calculations which include three- nucleon forces. 
Only such consistent calculations should be used to analyze the data in cases when both relativity 
and three-nucleon force effects are large. 

The paper is organized as follows. Sec. [TI] provides the conceptual basis for the choice of the 
momentum-space representation and the definition of spin in the relativistic context. In Sec. IIIII 
we summarize the formalism underlying relativistic three-nucleon Faddeev calculations with only 
nucleon-nucleon interactions, presented in detail in 
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121 • In Sec. lIVI we focus on the three-nucleon 



Faddeev equation with an included three-nucleon force and discuss two methods to compute matrix 
elements of the three-nucleon force in the partial wave basis used in our relativistic calculations. In 
Sec.|V]we apply our formulation to elastic nucleon-deuteron scattering and breakup and show and 
discuss the results. Sec. IVII contains our conclusions and summary. Appendixes lAl and [B] formulate 
three-nucleon forces in the momentum space representation adapted to Poincare invariance. 
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II. RELATIVISTIC DYNAMICS 



Relativistic invariance of a quantum theory means that the Poincare group (inhomogeneous 
Lorentz group) is a symmetry grou p of the theory. This requires the existence of a unitary repre- 



sentation of the Poincare group 



18|. The Poincare group has ten generators, six Lorentz generators 



J^'^, and four spacetime translation generators, P^. The dynamics of the system is given by the 
Hamiltonian, H = . The Lie algebra has two polynomial invariants. 



= -P^'P^ = W^Wf, (1) 

where is the Pauli-Lubanski vector {l^ 

:= -^e^^-f^yp^Jf,^. (2) 

It satisfies 

[p/^, W"] = 0, Pf'W^ = 0, [W, W] = ie^'^^^Pc.W/j. (3) 

Equation ([T|) implies that the Hamiltonian can be expressed in terms of the mass operator, H = 
VM2 + P2^ where .= p2^ Tj^^g^ 

given a representation for P, the dynamics is defined by the 
mass operator M, which plays the same role in Poincare invariant quantum mechanics as the center 
of mass Hamiltonian h = H — P"^ /2M does in Galilean invariant quantum mechanics. 

In the absence of interactions the mass operator M becomes the invariant mass operator Mq of 
three non-interacting relativistic particles. The full interaction is defined by 

V ■.= M- Mo. (4) 

For a system of three particles interacting with short-range interactions, two-body interactions 
are defined by 

^(ij)(k) '■= M(ij){k) - Mo (5) 

where M(jj)(fc) is obtained from M by turning off all interactions in M that involve particle k. The 
difference 

V4^:=V- V(i2)(3) - V(23)(l) - ^(31)(2) (6) 

defines a three-body interaction. With these definitions the mass operator has the form 

M = Mo + ^(12)0) + V^23)il) + ^(31)(2) + V4. (7) 



This has the same form as the non-relativistic three-body center of mass Hamiltonian with two 
and three-body forces, except the non-relativistic kinetic energy is replaced by the relativistic 
invariant mass of the non-interacting system. As in the non-relativistic case, bound and scattering 
eigenstates of this mass operator can be computed using the Faddeev equations with two and 
three-body interactions. For identical nucleons the coupled relativistic Faddeev equations can be 
replaced by a single equation. Details are discussed in the next section. 

In addition to the constraints imposed by discrete symmetries, translational invariance, and 
particle exchange symmetry, there are non-trivial constraints on the interactions due to both the 
Poincare symmetry and cluster properties. The constraints on the interaction due to Poincare 
invariance come from the commutator 

[P^,j''']=i6,kH, (8) 

which means that interactions appearing in H must be generated by the operators in the commu- 
tator. One way to satisfy the constraints due to Poincare invariance was suggested by Bakamjian 
and Thomas [201]. Their construction adds interactions to the mass Casimir operator that commute 
with the spin Casimir operator 

j2 := W^/M^. (9) 

The required interactions commute with and are independent of the total momentum and commute 
with the non-interacting three-body canonical spin operator. 

Spin is associated with rotational degrees of freedom that appear in the rest frame. Because 
the Lorentz boost generators, J*^*, do not form a closed sub-algebra, a sequence of Lorentz boosts 
that map the rest frame to the rest frame can generate a rotation. Thus in order to obtain a well- 
defined relativistic spin it is necessary to define a standard procedure for measuring the spin. This 
normally requires the specification of a special frame where spins can be compared (usually the rest 
frame) and a standard set of Lorentz transformations B^^{P)^y, parameterized by momentum, 
that transform arbitrary frames to the special frame. The three-body canonical spin is defined in 



terms of the Pauli-Lubanski vector by 



(OJc)^ := ^B-\Pr,W^ (10) 
where B~^{P)^y is the rotationless Lorentz transformation- valued function of the four momentum 

I 



B-\PT. :-- 



P^/M -P/M 



(11) 



This Lorentz transformation (jlip satisfies 

B-^PY^P" = {M, 0)'^. (12) 

Equations ([3]) and (jl2l) can be used to show that the components of jc satisfy SU (2) commutation 
relations. The spin (0,jc)^ is not a four vector because B~^{P)^y is a matrix of operators, rather 
than a constant Lorentz transformation. Under Lorentz transformation the canonical spin Wigner 
rotates 

(OJ^r :=ii^,(A,P)^(0,jer (13) 

where Rwci^, P)^u ■= {B^^{AP)ABc{P))^u- The spin Casimir operator = jc • jc is independent 
of the choice of boost (jlip used to define the spin. The non-interacting (kinematic) canonical spin, 
jco, is obtained from ^ by replacing M Mq, W -s- W^^ in ^ and M Mq in 1^. Thus, 
Poincare invariance can be satisfied provided the interactions y(ij){k) ^ commute with jco. 

The other non-trivial constraint on the interactions is imposed by cluster properties. The prob- 
lem arises due to the non-linear relation between the two-body interaction Vij in the two-body 
problem and the corresponding two-body interaction, in the three-body problem. Clus- 

ter properties relate ^(jj)(fc) to the Poincare generators for the interacting ij pair and spectator 
k. Unfortunately each 2-1-1 mass operator constructed by requiring cluster properties commutes 
with a different spin Casimir operator, which means that linear combinations of these interactions 
will break Poincare invariance. Coester [24] observed that these interactions could be replaced by 
phase-equivalent interactions that commute with jcO- These interactions are designed to satisfy 
cluster properties in the three-body rest frame. Using the Bakamjian-Thomas construction lin- 
ear combinations of the phase equivalent V(jj)(j;.)'s can be added in a manner that preserves the 
overall Poincare invariance. While these interactions do not lead to generators that satisfy cluster 
properties, cluster properties in the three-body rest frame and Poincare invariance of the S matrix 
ensures that the three-body S'-matrix retains cluster properties in all frames. 

To construct two-body interactions, the two-body interactions in the two-body problem that 
commute with the two-body canonical spin are replaced by phase equivalent two-body interactions 
in the three-body problem that commute with the three-body canonical spin. The phase equivalent 
interactions are identified in the rest frame of the three-body system. They are determined in all 
other frames by the requirement that the three-body spin remains kinematic (in the Bakamjian- 
Thomas construction this choice fixes the representation of the boost generators). 

To construct interactions that commute with the three-body kinematic canonical spin it is 
useful to introduce momenta and spin variables that have the same Wigner rotation properties 



as the three-body kinematic canonical spin. This is because the kinematic canonical spin can 
be constructed out of these degrees of freedom using conventional methods for adding angular 
momenta. 

The desired momentum operators are the relativistic analog of Jacobi momenta. In the non- 
relativistic case Jacobi momenta can be defined using Galilean boosts to the two and three-body 
rest frames. In the relativistic case the Galilean boosts are replaced by the rotationless boost pip 



and the relevant Jacobi momenta are 



21 



q1 = B-\P)\p'; P''=pUpf^+P^ qt = i^qf+m'^,q^ (14) 
k^, = B-\q,,nq'^ q^:=q^ + q^. (15) 



In terms of these variables 
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Mo = J2^m^+qf = Jml^ + qi + ^rn?+ql (16) 



i=l 

where the two-body invariant mass is 



"^uo = \J-{qt + qjY (g» + )/. = ^m^^ + + + . (17) 
The vector variables satisfy 

3 

^ q, = k,j + kji = 0. (18) 

i=l 

The relevant property of these momentum vectors is that they experience the same Wigner rotations 
as the three-body kinematic canonical spin ()13p . 

q1 ^ qf = {B-\kP)KB,{P)B-\P)Y,p\ = R^,{A, Py.q^ . (19) 



Similarly, 



k1^ ^ k>t!j = B-\q[^Y,qf = {B-\RUK,P)q,,)R^,{k,P)Y,q'i 



{B-\R^,{A, P)qij)R^,{k, P)B,{qij)YuKj = R^c{K PYi^Kj (20) 



where the last line follows from the property of the rotationless boosts (jlll) that the Wigner rotation 
of a rotation is the rotation j2l| 

R^Mpru = R^. (21) 
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for any P. Thus the and kjj all undergo the same Wigner rotations as the three-body kinematic 
canonical spin. 

Next we introduce single-particle spins with the same property. Single-particle canonical spins 
can be constructed from single-particle Poincare generators using 

{o,j,,r ■■= ^B^Hp^r.wr (22) 

where the operators on the right side of (j22p are constructed by replacing all of the three-body 
generators in ([2]), ()10p and (llip by the corresponding one-body generators. 

Under kinematic Lorentz transformations the single-particle canonical spins experience Wigner 
rotations, Rwc{^,Pi)^u, that depend on the single-particle momenta. These rotations differ from 
the Wi gner rotations experienced by cjj, kjj and jco- This can be changed by introducing new 
single-particle spin operators that replace the rotationless boost in (p2]) by a two step boost, 

B-\p,r, ^ {B-\q,)B^\P)r,. (23) 

These two boosts agree when P = 0. Note that both of these boosts transform — t- (m, 0,0,0)^, 
so they differ by momentum dependent rotations. We call these spins three-body constituent spins 
to distinguish them from single-particle canonical spins. The constituent spin operators are defined 



by 



22 



(0,j3c.^)'^ := -{B~\q,)B-\P)Y,W^. (24) 
m 



When P = 0, B^^{pi) B^^{qi) which means that single-particle canonical spins and three-body 
constituent spins agree in the three-body rest frame. For a three-body system the total spin is 
identified with total angular momentum in the three-body rest frame, which is the sum of the 
single-particle angular momenta. The angular momentum of a single particle in the three-body 
rest frame is the sum of the single-particle constituent spin and a contribution from the single 
particle orbital angular momenta. 

A calculation, using the property (j2ip . shows that under Lorentz transformations 



(0, j'3,,J'^ := i?^e(A, Pru{0,hcsir , (25) 

Wigner rotates with the same rotation as the vectors and k,j and the three-body kinematic 
canonical spin. The three-body kinematic canonical spin is the sum of the orbital angular momenta 
associated with q^. and kjj and the single-particle three-body constituent spins. The requirement 
that an interaction commutes with the kinematic three-body canonical spin is equivalent to the 



requirement that the interaction have a rotationally invariant kernel when expressed in terms of 
these variables. Thus the required interactions in the Bakamjian-Thomas construction are given 
by kernels of the form 

(P,qi,kjfc,//i,/i2,A^3|^|P',qi,kjfc;/^'l>/^2>A*3) = ^(P - P')(qi;kjfc> W,M2,M3||^||qi,kjfe,At'l,At2,/^3) 

(26) 

where the reduced kernel is a rotationally- invariant function of qi,kjfc and the three-body con- 
stituent spins. 

Two-body interactions in the two-body problem vi2 have a similar form 

(Pl2,ki2,/il,^2bl2|Pi2>ki2,/i'i,^2> = '^(^12 " P12) (^12 , W , /^2 1^12 ||ki2 , /u'l , /is) (27) 

where 

k'^^ = B-\p,+p2ruPf (28) 

is the two-body relative momentum and the magnetic quantum numbers are associated with the 
two-body constituent spins 

{0,hcsir = ^{B-\kij)B-Hpi+Pj)ruW,r (29) 

When these interactions are embedded in the three-body Hilbert space the kernels (1271) are 
replaced by kernels that are rotationally invariant functions of the three-body Jacobi momenta 
and the three-body constituent spins. In order to satisfy cluster properties kjj given by ()28p is 
replaced by the kjj given by ([15]) , the pi are replaced by the corresponding Qi , and the two-body 
constituent spins ([29]) are replaced by 

(0^2(3)^^)^^ = UB-Hhj)B-\q, + q,)B-\P)r,Wr. (30) 

These operators represent two-body constituent spins in the three-body rest frame. They agree 
with the two-body constituent spins (j29p that they replace in the three-body rest frame, but are 
defined so they remain unchanged by canonical boosts out of the three-body rest frame. This 
ensures that they undergo the same Wigner rotations as the kinematic three-body canonical spin 
under kinematic Lorentz transformations. Thus, the kernels (|27p are related by 

(P, Qj, kjfc, fJ.l,p,2,tJ'3\Vjk\'P', <ii, kj-fc, Pi,fi2,fJ-3) = 

5(P - P')<5(q. - c{i)S^^^, J2 [Bc{Qj)Bc{qj + 9fc)i?c(%fc)]I?i{L [BMB^Qj + qk)Bci-k,k)] x 
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Dl/^^,\B-\^,)B-\q'^+q',)B-\0^^^ (31) 

Here the unbarred magnetic quantum numbers are three-body constituent spins while the barred 
magnetic quantum numbers are the two-body constituent spins in the three-body rest frame. 

Even though the spins in (|30|) transform the same way as the three-body constituent spins, they 
differ from the three-body constituent spins (j24p by the Wigner rotation 

{0,h(3)icsr = {{B-\ki,)B-\qi + qi)B-\qi)Y,{Q,hicsY- (32) 

When the two-body interactions are embedded in the three-body system the spins are identified 
with the two-body constituent spins in the three-body rest frame, as would be expected by cluster 
properties, but in other frames they are defined to remain unchanged with respect to canonical 
boosts out of the three-body rest frame. The Wigner rotations (|44p and (|A2p arise because the 
two-body subsystem is moving in the three-body rest frame; however because the Wigner rotations 
in (j32p are functions of the qi rather than the pj, both spins in (j32p undergo the same Wigner 
rotations under kinematic Lorentz transformation. Because of this it is also possible to construct the 
three-body canonical spin using partial wave methods directly in a mixed representation involving 
the barred spins in the interacting pair and the unbarred spin for the spectator. In the mixed 
representation the two-body interaction in the three-body Hilbert space has the simple form 

(P, qi, kjfc, /i2, AisbjfclP', Qi, k^-fc, f^'i, (^'2, f^'s) = 

5{P - P')5{cii - c^^)5f,^^>,(kJk,^ij,^lk\\vjk\\^k^P''j^P''k) (33) 

For the two-body problem in the three-body Hilbert space it is advantageous to use ([33]) because 
spins (|3L)p do not require Wigner rotations. However, with this choice each interacting pair of 
particle must be treated using a permuted basis which requires Wigner rotations in the permutation 
operators. The three-body forces are naturally expressed by a rotationally invariant kernel in the 
three-body constituent spins. When they are transformed to a mixed basis that involves the spin 
(|30p for one pair, then it is necessary to transform two of the three-body constituent spins with the 
Winger rotations in (j32p . The calculations performed in this work use a partial wave projection 
of the mixed basis (j33p . although the Wigner rotations in the three body-interaction are not yet 
included. 
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III. RELATIVISTIC THREE-NUCLEON FADDEEV EQUATIONS WITH NUCLEON- 
NUCLEON FORCES 



The nucleon-deuteron scattering with neutron and protons interacting through only a nucleon- 
nucleon interaction vmmis described in terms of a breakup operator T satisfying the Faddeev-type 



integral equation 



23 



2M 



T\cP) = tP\cl))+tPGoT\<p). (34) 

The two-nucleon t-matrix t is the solution of the Lippmann-Schwinger equation with the interaction 
vnn. The permutation operator P = P12P23 + P13P23 is given in terms of the transposition 
operators, Pij, which interchanges nucleons i and j. The incoming state \(j)) = |qo)|i?^d) describes 
the free nucleon-deuteron motion with relative momentum qo and the deuteron state vector \(j)d)- 
Finally Go is resolvent of the three-body center of mass kinetic energy. Transition operators for 
the elastic nd scattering, U, and breakup, Uq, are given in terms of T by [2J 

U = PGq^ + PT , 



f/o = (1 + P)T . (35) 

This is our standard nonrelativistic formulation, which is equivalent to the nonrelativistic three- 
nucleon Schrodinger equation plus boundary conditions. The formal structure of these e quat ions 
in the relativistic case remains the same but the ingredients change. As explained in [25] the 
relativistic three-nucleon rest Hamiltonian (mass operator) has the same form as the nonrelativistic 
one, only the momentum dependence of the kinetic energy and the relation of the pair interactions 
in the three-body problem to the pair interactions in the two-body problem change. Consequently 
all the formal steps leading to (I34p and (|35p remain the same. 

The free relativistic invariant mass of three identical nucleons of mass m has the form [1^ (see 
Eq.dlSI)) 



Mo = ^ml^Q + q^ + ^m'^ + q'^ (36) 

with spectator momentum q := qi and the free two-body mass operator m23o expressed in terms 
of the relative momentum k := in the 2 — 3 center of momentum frame by (see Eq. ()17p ) 



771-230 = 2VW+rn? = 2um{k) ■ (37) 

As introduced in [22] and in Eq.Q the pair forces in the relativistic three-nucleon 2-1-1 mass 
operator are related to the two-body forces in the two-body problem, Vij, by 



= V ^^ijo + Vij)"^ + 9^ - V "ifjo + 5^ ' (38) 
11 



where V = V{q'^) reduces to the interaction v for q = 0, which acts in the two-body center of 
momentum frame. The momentum dependence ensures that the resulting three-nucleon scattering 
matrix satisfies space-hke cluster properties in all frames [22l |. 



The transition matrix t that appears in the kernel of the Faddeev equation (|34p is obtained by 
solving the relativistic Lippmann-Schwinger equation as a function of 



^{2u;m{k')r + q'- ^{2u:^{k")f + + ie 
The input two-body interaction V is computed by solving the nonlinear equation 13 1 



{ V"^?jO + %)(fc)} + %)(fc) = "imvNN , (40) 

where v^n is a nonrelativistic nucleon-nucleon potential fitted to the nucleon-nucleon data basis 
and where anticommutator {A, B} = AB + BA. In case of q = that equation reduces to a 
nonlinear equation for the relativistic two-body interaction v. Therefore the problem of refitting 
all two-nucleon data when changing from a nonrelativistic to a relativistic Lippmann-Schwinger 
equation is avoided. The nonlinear equation (j40p can be solved by iteration 13] • An alternative 
approach to determine t(k, k';g^) is described in |l4|. 

The new relativistic ingredients in ()34p and (j35p will therefore be the i-operator (I39p (expressed 
in partial waves) and the resolvent of the three-nucleon invariant mass 

with Mq given by (j36p . E is the total three-nucleon invariant mass expressed in terms of the initial 
neutron momentum qo relative to the deuteron by 



E = ^Ml+q,^ + ^m^ + q^\ (42) 

with Md the deuteron rest mass. Related to the choice of the permutation operator P the pair i — j 
is chosen as 2 — 3. 

Currently the Faddeev equation (j34p in its nonrelativistic form is numerically solved for any 
nucleon-nucleon interaction using a momentum space partial-wave decomposition. Details are 
presented in [2^3]. Projecting on such a basis turns it into a coupled set of two-dimensional 
integral equations. As shown in [ll|, [l^, in the relativistic case we can keep the same formal 
structure, though the permutation operators are replaced by the corresponding Racah coefficients 
for the Poincare group which include both Jacobians and Wigner rotations that do not appear in 
the nonrelativistic permutation operators 



24 
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In the nonrelativistic case the partial- wave projected momentum-space basis is 



\pq{ls)j{X^)IJ{t^)T), 



(43) 



where p and q are the magnitudes of standard Jacobi momenta (see 2J, |26|]), obtained by trans- 
forming single particle momenta to the rest frame of a two- or three-body system using a Galilean 
boost, and {ls)j are two-body quantum numbers with obvious meaning, (A1/2)I refer to the third, 
spectator nucleon, taken as the nucleon 1 and described by the momentum q, J is the total three- 
nucleon angular momentum and the rest are isospin quantum numbers. In the relativistic case this 
basis is replaced by the Poincare irreducible states defined as |]J 



(Pl,A^'l,P2,/^2,P3,/^3l('^,9)P = 0,^; A,I, j23,A:23,^23,'S23) = 

(5 - qi - q2 - qsj^TTT r 2 U 

(2' -,H3\s,Hs)il,fJ'hS,Hs, |j,/ij)(A,/iA, -,n'i\I,fIl){j,lij,I,Hl\J,tl) 

MSA's 

Yx^Ml)Yl^.AH^2.m))D\^,[Ru,c{B,{-^^^ (44) 

where N(q2, qs) is given by (|A3p in Appendix [Aland k{q2, qs) by (|Bip in Appendix[Bj These states 
are labeled by the same quantum numbers as the corresponding non-relativistic basis states. 

The basis states (I44p are used for the evaluation of the partial wave representation of the 
permutation operator P with Wigner rotations of spin states for nucleons 2 and 3 included. In the 
relativistic case we adopt the following short-hand notation for the Poincare irreducible three-body 
states, which also includes isospin quantum numbers coupled in the same order: 

\k, q, a) := \kqils)j{\, ^)IJ{t^)T) = |(J, q)P = 0, fi; A, /,i23, ^23, ^23, ^23)1(4)^) • (45) 
Equipped with that, projec ting onto the basis states \k,q,a) one encounters, using the non- 



relativistic notation of Ref. 



26| 



i{kqa\P\k'q'a)i 



i{kqa\k' q'a)2 +1 {kqa\k' q a)^ = 2 i{kqa\k' q a 



(46) 



This is evaluated by inserting the complete basis of single-particle states |pi, P2, /^2, P3i /^s) and 
using (j44p . It can be expressed in a form which resembles closely the corresponding non-relativistic 



expression [24 



i{k q a\ P \k' q d)\ 



1 b{k- vri) b{k' - 7r2) 
ax- 



-1 



A;2 



iVi(Q,Q',x) iV2(g,g',x) 



Gaa'iq^ q' 



(47) 
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where all ingredients are defined in Appendix IA2I of Ref . [10] . The rotational invariance of the 
nucleon-nucleon interaction in this basis ensures that all three nucleon-nucleon interactions com- 
mute with the spin Casimir operator of the non-interacting three-nucleon system. This allows the 
interactions to be added in a manner that preserves the underlying Poincare symmetry. 

Due to the short-range nature of the nucleon-nucleon interaction it can be considered negligible 
beyond a certain value jmax of the total angular momentum in the two-nucleon subsystem. Gen- 
erally with increasing energy jmax will also increase. For j > jmax we set the t-matrix to zero, 
which yields a finite number of coupled channels for each total angular momentum J and total 
parity vr = (— )'^'^ of the three-nucleon system. To achieve converged results at incoming nucleon 
laboratory energies below ~ 250 MeV all partial wave states with total angular momenta of the 
two-nucleon subsystem up to jmax = 5 and all total angular momenta of the three-nucleon system 
up to J = 25/2 must be taken into account. This leads to a system of up to 143 coupled integral 
equations in two continuous variables for a given J and parity. For the details of the numerical 
performance we refer to 



11 



24 



26( 1 . The solution of these equations can be used to construct an 



exactly Poincare invariant scattering operator. 



IV. RELATIVISTIC THREE NUCLEON FADDEEV EQUATIONS WITH A THREE- 
NUCLEON FORCE INCLUDED 

In the standard nonrelativistic formulation when in addition to pairwise interactions i^attv be- 
tween three nucleons also a three-nucleon force is included, a new term V4 appears in a potential 
energy of the three-nucleon system 

= + V}^^ + V}^^ . (48) 

Each v}^^ is symmetric under exchange of the nucleons j and k {i,j,k = 1,2,3 and j ^ i ^ k). 
In the 27r-exchange three-nucleon force V^^^^ is a contribution to the three-nucleon potential from 
(off-shell) rescattering of a pion on nucleon 1. 

When a three-nucleon force is acting then on top of rescatterings among three nucleons induced 
by pairwise forces only, which are summed up in integral equation (|34p . additional rescatterings 
induced by three-nucleon force and nucleon-nucleon force appear. 

Therefore Faddeev equation ()34p changes to 

T|0)= tP\(j)) + {l + tGo)vi'\l + P)\^) +tPGoT\(l)) + {l + tGo)vi'^ {l + P)GoT\^) , (49) 
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with one new contribution in the leading term and in the kernel [2J, |27|] . While the breakup transi- 
tion operator Uo preserves its form (|35p . in the elastic scattering operator U two new contributions 
appear 



24 



U = PGq^ + Vi^^ (1 + -P) + PT + Fi'^ (1 + P)GoT 



(1) 



(50) 



The second term is due to a single interaction of three-nucleons via a three-nucleon force and the 
fourth results from rescattering among three nucleons induced by two- and three-nucleon forces 
with a three-nucleon force as the final interaction. 

After projecting on a partial-wave momentum-space basis equation ()49p becomes a system of 
2-dimensional coupled integral equations which can be solved numerically exactly for any nuclear 
force. Since the three-nucleon force is short-ranged its inclusion needs to be carried through only 
for all total angular momenta of the three nucleon system up to J = 13/2. As mentioned in section 
mil the longer ranged two- nucleon interactions require states up to J = 25/2. For details of the 
formalism and numerical performance in case of the nonrelativistic formulation when three-nucleon 
force is included we refer to Refs. 



231, 



24 



28| 



For relativistic calculations without a three-nucleon force, briefly described in previous section, 
the details of the numerical treatment are given in 
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121. When a three-nucleon force is added 



two new terms in ()49p contain the free three-nucleon propagator Gq. Since in the basis |k, q > (see 
Appendix [B]) the three-nucleon invariant mass Mq is diagonal, Gq is given by 

1 



(k,q|Go|k',q') = 5(k-k05(q-q') 



(51) 



E - y^m'^ + - ^J'i{k'^ + m^) + + ie 
That means that performing integrations over momenta k' and q' in the intermediate states 
\k' , q' , a') during the calculation of matrix elements for these new terms, the simple pole singularity 
occurs for momenta q' < qmax at k' = k^, where qmax is given by the total three-nucleon center of 



momentum energy E through E = + q^ax + a/"^^ + Qmax • 

For a given g'-value the momentum k^ is the solution oi E = J 4:{m? + /cq) -\- q'"^ + y m? + q^ 



v2 



The treatment of that pole, as well as of the deuteron bound state pole in T, which occurs at 
= qo for channels a' containing the deuteron quantum numbers, was done using subtraction 



method 

The nonrelativistic treatment of (j49p requires matrix elements of V^\l + P) calculated in a 
partial-wave basis with standard Jacobi momenta: (p, q, a\V^'^\l + P)\p' , q' , a'). In the relativistic 
calculations, however, one needs them in the new, relativistic basis \k,q,a). The generation of 
three-nucleon force partial- wave matrix elements is the most time consuming part of three-nucleon 
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continuum Faddeev calculations. One way to reduce the computer time is to perform a trans- 
formation of the existing, standard Jacobi momenta matrix elements to the relativistic basis. In 
Appendix |X] we give the expression ()A17p for such transformation which is valid in the general 
case, when in addition to boost also Wigner spin rotations are taken into account. The complex 
structure of that transformation, where in addition to the summation over numerous intermedi- 
ate states with geometrical coefficients, also involved are two integrations and two interpolations 
over the momenta p and p' , prevents, due to the large amount of computing time and computer 
resources required, the application of that transformation in fully converged calculations. 

It seems thus unavoidable that in order to get matrix elements {k,q,a\V^^\k\q\a') one must 
start from a commonly given expression for a three-nucleon force in terms of individual nucleons 
momenta and their spin and isospin operators an d to apply to that expression the recently proposed 
automatized partial wave decomposition 



291, 



SOj. To that aim we derived in Appendix [B] relation 



()B4p which allows to express matrix element of a three-nucleon force in a 3-dimensional relativistic 
basis (k, q|V^'""'^^|k', q') by its matrix element in the individual nucleons momenta basis. In this basis 
q and k undergo identical Wigner rotations under kinematic boost of the three-nucleon system. The 
nucleon spins are defined to be the three-body constituent spins (canonical spins measured by using 
a rotationless boost to the three-body center of momentum frame). To use them in the \k,q,a) 
basis the spins for the a pair must be Wigner rotated before they are coupled. The alternative is 
to use the representation where three spins are three-body constituent spins; in this case all three 
of the two-body interactions will have Wigner rotations that convert the two-body constituent 
spins in the three body rest frame to three-body constituent spins. The three-nucleon force will 
have no Wigner rotations. In this representation all of the spins can be coupled using standard 
partial wave methods. For our calculations we work in the \k,q,a) basis, but do not account 
for the Wigner rotations in the three-nucleon interaction for the reasons discussed in the previous 
paragraph. This allows us to treat the spins in the three-nucleon force using conventional methods. 
This assumes the Wigner rotations in the three-nucleon force can be neglected. Neglecting these 
Wigner rotations has no effect on the relativistic invariance or S-matrix cluster properties. 



V. RESULTS 



To study the importance of a consistent treatment of both relativity and a three-nucleon force 
we numerically solved the three-nucleon Faddeev equations for neutron-deuteron scattering at the 
neutron laboratory energies E^^^ = 70, 135, 200 and 250 MeV. As dynamical input we took the 



16 







31 



33]. The 



nonrelativistic nucleon-nucleon potential CD Bonn [4] and TM99 three-nucleon force 
cut-off parameter A of that three-nucleon force was adjusted to A = 4.469 in units of the pion mass, 
mjr, to give, together with the CD Bonn potential, the experimental binding energy of '^H. At each 
energy we generated solutions of nonrelativistic and relativistic three-nucleon Faddeev equation, 
without and with TM99 three-nucleon force included. For relativistic case we produced, starting 
from the CD Bonn potential and solving nonlinear equation (j40p at the required spectator nucleon 
momenta q, the relativistic, on-shell equivalent interaction with boost effects incorporated exactly. 
That interaction served as dynamical input to calculate, using the relativistic Lippmann-Schwinger 
equation (j39l) the relativistic off-shell t-matrix t that appears in Faddeev equations. 

Since in [l2] it was found that effects of Wigner spin rotations are practically negligible in 
the studied energy range, we neglected them in the present study. When performing relativistic 
calculations with three-nucleon force included one requires matrix elements of the TM99 three- 
nucleon force in a relativistic momentum space basis, where the relative momentum of two nucleons 
in their cm. system, k, replaces standard Jacobi momentum p. That momentum k together with 
spectator nucleon momentum q, equal in magnitude and opposite to the total momentum of the 
free pair in three-nucleon center of momentum system, unambiguously define the configuration of 
three nucleons. Since it is the region of small and not large momenta which is most important when 
solving Faddeev equations it seems reasonable to assume that the momenta k and p do not differ 
substantially. Therefore, in order to avoid calculations of the TM99 three-nucleon force matrix 
elements in a relativistic basis \k,q,a) we assumed, that the matrix elements in a relativistic and 
nonrelativistic bases are equal: 

{k,q,a\vj^^\k' ,q ,a) = (p = k,q,a\v}^^\p = k' ,q ,a) . (52) 

That assumption allowed us to use the existing matrix elements of the TM99 three-nucleon force. 

To check quality of the approximation ()52p we compared the matrix element of the TM99 
3NF in the relativistic basis, (fc, q, a\v}^^\k' , q', a'), calculated according to (|B4|1 and using au- 



tomatized partial wave expansion of Ref. 30| (what corresponds to neglection of Wigner spin 
rotations in ()B5P ). with the corresponding matrix element in the standard, nonrelativistic basis, 
{p,q,a\v}^^\p',q',a'), at a number of the spectator momentum values. In Figs. [1] and [2] we ex- 
emplify the typical behavior showing at a number of q' values and at a fixed p = k, taking two 
different values of q, the k'(= p') dependence of these matrix elements for a particular channel 
a = a' = 1(00)0(0^)^(1^)^ >. As expected, clear differences between these matrix elements occur 
only at very large values of the spectator momentum q, where magnitudes of these matrix elements 
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are small. This justifies application of the approximation (j52p in the present study. 

The approximation (|5'2p can be investigated also directly for the three-dimensional matrix ele- 
ments, comparing (k, q | V^^^^ |k',q') and (p,q | V^^"^^ |p',q')- The connection between these matrix 
elements is given by ()B4p in Appendix [Bl They depend on momentum vectors and spin-isospin 
quantum numbers in the initial and final state. In Fig. Owe show a particularly simple case, where 
t = t' = 0, all four momenta are parallel to the unit vector ^--i=, -i=, -i= j and all spin magnetic 
quantum numbers are equal ^. We display (k, q | vj^^ |k, q) and (p, q | v}^^ |p, q) for several q val- 
ues as a function of k. We see how the difference develops gradually with increasing q, resembling 
the picture seen for partial wave decomposed matrix elements. 

Transition amplitudes for elastic neutron-deuteron scattering and breakup based on that set of 
four solutions of three- nucleon Faddeev equations, are used to predict numerous observables for 
both reactions. By comparing these observables conclusions on how strongly three- nucleon force 
effects depend on relativity were drawn. In the following subsections we show and discuss results 
for the cross section and numerous spin observables, separately for elastic scattering and breakup 
reactions. 



A. Elastic scattering 

At higher energies of the incoming nucleon three-nucleon forces play significant role in deter- 
mining the angular distribution of the elastic neutron deuteron scattering. The clear evidence 
of three-body force effects start to develop at El^'^ ~ 65 MeV for scattering angles close to a 
minimum of the cross section, which at 65 MeV occurs at Oc.m. ~ 105° [g, J With increasing 
energy of the three-nucleon system not only the magnitude of predicted three-nucleon force effect 



uences the cross section in a wider range of angles, which at 250 MeV covers 
The standard 27r-exchange three-nucleon forces, such as TM99 33( or 



increases but it also inf 
90" < Ocm. < 180° 

Urbana IX [s^, are able to account for existing discrepancies between theoretical cross sections 
obtained with realistic nucleon-nucleon potentials and data only up to El^^ ~ 135 MeV. Data at 
larger energies in a region of angles ranging from the cross section minimum up to 180° are drasti- 
cally underestimated even when 27r-exchange three-nucleon forces are included in the calculations. 
This is exemplified on Fig. [H where solid (red) lines are nonrelativistic predictions based on the 
CD Bonn potential alone and dotted (blue) lines are results obtained when the CD Bonn potential 
was combined with the TM99 three-nucleon force. 

Since effects of relativity for predictions based on two-nucleon forces only are restricted to very 
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backward angles 6c.m. > 160" [111 (see also Fig. d] where dashed (blue) lines are relativistic predic- 
tions based on the CD Bonn potential), the drastic discrepancy between data and theory seen at 
250 MeV would indicate that at such large energies shorter-ranged three-nucleon force components, 
not taken into account in these calculations, start to play significant role. The possibility, that in- 
cluding such three-nucleon force contributions would indeed help to improve description of the cross 
section data is further supported by an interesting pattern revealed when the TM99 three-nucleon 
force is included into relativistic calculations. Namely, when a consistent treatment of relativity 
and a three-nucleon force as described in the present study is made, then the resulting changes of 
the cross section are not a simple incoherent sum of effects due to relativity, seen when two-nucleon 
forces alone are acting, and three-nucleon force effects found in nonrelativistic calculations. The 
relativity modulates effects exerted by the TM99 three-nucleon force on the cross section found 
in nonrelativistic calculations and the magnitude of this modulation depends from the scattering 
angle. While at backward angles the nonrelativistic cross section with a three-nucleon force in- 
cluded is further enhanced by relativity, in a region of center of momentum angles near the cross 
section minimum the magnitude of three-nucleon force effects seen in nonrelativistic calculations 
is strongly reduced by relativity (dashed-dotted (brown) lines in Fig. |3| . 

Also elastic scattering polarization observables reveal such incoherent and angle-dependent mod- 
ulation of three-nucleon force effects by relativity. The details, however, depend on the particular 
spin observable under study and every conceivable scenario can be found. 

For elastic scattering spin observables effects of relativity, when only two nucleon forces are 
acting, were found to be small 1J|. It is exemplified by nearly overlapping solid (red) and dashed 
(blue) lines in Figs. [5)1151 Adding three-nucleon force in nonrelativistic calculations leads to sub- 
stantial effects for some polarization observables, especially at higher energies 0)0] • The resulting 
picture, however, is quite complex. Some of those three-nucleon force effects are supported by the 
data. For some observables they deteriorate the data description. 

For tensor analyzing powers Axx, ^yy and A^z relativistic effects are non- negligible even at 
70 MeV (see Fig. [5]) and clearly increase with increasing energy as seen in Figs. [6l [7] and [9j When 
three-nucleon force is added in the relativistic calculations the resulting effect depends on the 
observable and the energy. 

For Axz large three-nucleon force effects remain. At 70 MeV and 135 MeV they are prac- 
tically identical in magnitude to three-nucleon force effects found in nonrelativistic calculations 
and nonrelativistic and relativistic predictions for A^z at these energies are practically overlapping 
(see dotted (blue) and dashed-dotted (brown) lines in Figs. [5l [6] and [9]). At 200 MeV, however. 
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adding three-nucleon force in relativistic calculations leads to angle dependent modulations of the 
magnitude of three-nucleon force effects, similar to that found for the cross section (see Fig. E]). 

For Axx a drastically different scenario occurs. Large effects of the TM99 three-nucleon force 
are seen for that observable in nonrelativistic calculations at 70 MeV and 135 MeV in wide range 
of angles and they practically vanish when relativity is included. As a result the dashed-dotted 
(brown) line practically overlaps with pure two-nucleon relativistic and nonrelativistic predictions 
(see Fig. [5] and [6]). 

For Ayy (Fig. [5] and [6|) the large effects of the three-nucleon force seen in nonrelativistic calcu- 
lations are simply reduced by relativity. For the tensor analyzing power A^z, for which data exist 
only at 135 and 200 MeV, the influence of relativity induces both modulation and reduction of 
nonrelativistic three-nucleon force effects (Fig. [7]) . 

The TM99 three-nucleon force acts differently on the nucleon, Ay{N), and deuteron, Ay(d), 
vector analyzing powers. While three-nucleon force effects for Ay{N) are rather small even at 
250 MeV (Fig. El and 15]), for Ay{d) they are significant (Fig. [5l [6] and [8]). For Ay{N) and Ay{d), 
but more clearly displayed due to larger effects for the deuteron vector analyzing power, both 
reduction and modulation of nonrelativistic three-nucleon force effects by relativity was found. 
That reduction and modulation depend on angle and energy. 

A similar picture was found for numerous spin correlation coefficients, as exemplified by different 
theoretical predictions shown in Figs. [9]|13[ Again all scenarios are available: total reduction by 
relativity of large three-nucleon force effects seen in nonrelativistic calculations (e.g. Cx,x at 135 
and 200 MeV for 120° < 9c.m. < 150° in Fig. [JOl Cy^y at 135 and 200 MeV at 120° < 9c.m. < 
150° in Fig. [TTI) . practically the same three-nucleon force effects in nonrelativistic and relativistic 
calculations {Cz,z at 135 MeV in Fig. [9l Cxz,y at 135 and 200 MeV in Fig. [12]), angle dependent 
modulation of nonrelativistic three-nucleon force effects by relativity {Cx^z at 135 MeV in Fig. 
Cz,x at 135 and 200 MeV in Fig. lOl Cxy,x and Cyz,x at 135 and 200 MeV in Fig. \T3\i. 

The polarization transfer coefficients are not exceptions; also for them a similar complex influ- 
ence of relativity on nonrelativistic three-nucleon force effects have been found as shown in Figs. [T^ 
andll5[ 

The comparison of nonrelativistic predictions based on 27r-exchange three-nucleon force's re- 
vealed for spin observables a complex, angle and energy dependent pattern of discrepancies between 
data and theory The nontrivial interplay between the 27r-exchange three-nucleon 



forces and relativity suggests that the inclusion of further three-nucleon force mechanisms, like 
forces of shorter range, is needed to improve the description of elastic scattering polarisation data. 
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B. Breakup 



Theoretical study of exclusive breakup reaction performed at different incoming nucleon energies 
revealed regions of breakup phase-space where large three-nucleon force effects have been found 



38l |. The effects, similarly to elastic scattering, generally increase with energy. With increasing 



energy also the effects of relativity increase 
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16[, revealing for exclusive breakup cross section 



a characteristic pattern when viewed as a function of the angles of detected nucleons. Largest 
effects where found when two of three outgoing nucleons are detected coplanarly on both sides of 
the beam. Keeping one of the detectors at a constant position and changing the polar angle of 
the second, regions of phase space were found in which nonrelativistic breakup cross section was 

n 

increased or decreased by relativity [16]. In these specific configurations effects of three-nucleon 
force's on breakup cross section, both in norelativistic as well as in relativistic calculations, are 
practically negligible (see Fig. [TBI) . 

Due to richness of the breakup phase-space also geometrical configurations can be found where 
both, three-nucleon force and relativistic effects are significant. Exclusive cross sections in some of 
these configurations are shown as a function of the laboratory energy of one of the outgoing and 
detected nucleons in Fig. [T7] for neutron-deuteron breakup at 200 MeV. It is seen that including 
relativity reduces slightly the magnitude of three-nucleon force effects observed in nonrelativistic 
calculations. 

Relativity changes also the magnitude of three-nucleon force effects seen in nonrelativistic calcu- 
lations for breakup polarization observables. We exemplify that in Fig.llSlatjthree configurations of 
exclusive dp breakup at E"^"^ = 270 MeV, for which data have been taken [17]. Again, influence of 
relativity on magnitude of three-nucleon force effects change with configuration as shown in Fig. [18] 
along the S-curve arc length. Especially interesting is the case of polarization-transfer coefficient 
from the deuteron to the nucleon, Kl^'y, for which inclusion of TM99 three-nucleon force changes 
completely the S-dependence found in case when only two nucleon-forces were acting. The effect 
of three-nucleon force is further modified slightly by relativity resulting in a better reproduction of 
data. 

VI. SUMMARY AND OUTLOOK 

We extended our relativistic formulation of three-nucleon Faddeev equations to include also 
three-nucleon force. The relativistic features are the relativistic form of the free propagator, the 
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change of the nucleon-nucleon potential caused by the boost of the two nucleon subsystem, and the 
modification of the permutation operators. In present study we neglected Wigner spin rotations 
induced by these boosts since investigations based on two-nucleon forces only have shown that 
their effects are negligible. For the momentum-space basis we used the relative momentum of two 
free nucleons in their cm. system together with their total momentum in the three nucleon cm. 
system, which in this frame is the negative momentum of the spectator nucleon. Such a choice 
of momenta is adequate for relativistic kinematics and allows to generalize the nonrelativistic 
approach used to solve the nonrelativistic three nucleon Faddeev equation to the relativistic case 
in a more or less straightforward manner. That relative momentum in the two-nucleon subsystem is 
a generalization of the standard nonrelativistic Jacobi momentum p. We numerically solved these 
equations for neutron-deuteron scattering including relativistic features and/or three-nucleon force 
at the neutron lab energies El^^ = 70, 135, 200 and 250 MeV. As dynamical input we took the 
nonrelativistic nucleon-nucleon potential CD Bonn and generated in the two nucleon center of 
momentum system an exactly on-shell equivalent relativistic interaction. As a three-nucleon force 
we took the 27r-exchange TM99 force. 

By comparing our relativistic calculations without and with the three-nucleon force included 
we studied influence of relativity on three-nucleon force effects. In studies with two-nucleon forces 
only it was found that significant relativistic effects for the elastic scattering cross section appear at 
higher energies and they are restricted only to the very backward angles where relativity increases 
the nonrelativistic cross section. At other angles the effects are small. Also for spin observables, 
analyzing powers, spin correlation coefficients and spin transfer coefficients, no significant changes 
due to relativity have been found when only two-nucleon forces were acting. The similar picture 
was found for breakup, however, in that case significantly larger effects for the cross section in 
specific regions of the breakup phase-space have been found. 

The results obtained in the present study document that this picture changes dramatically when 
in addition to the two-nucleon force in a relativistic treatment also a three-nucleon force is acting. 
For the elastic scattering large changes of the cross section at higher energies, caused by three- 
nucleon force in large region of angles ranging from around minimum of the cross section up to very 
backward angles, arc further significantly modulated by relativity. Also such modulation in a large, 
similar to that for the cross section, range of angles have been found for numerous polarization 
observables. In that case every conceivable scenario of modulations was observed: from wiping 
out large three-nucleon force effects found in nonrelativistic calculations to their modulations with 
energy and angle, with strong amplification or reduction of their magnitude. Also for exclusive 
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breakup cross section and polarisation breakup observables in some geometries the relativity influ- 
ences effects induced by three- nucleon forces. Thus also for that reaction the relativistic treatment 
when three-nucleon forces are acting is required for proper interpretation of data. 

The comparison of our nonrelativistic theory with existing elastic scattering cross section and 
polarisation data exhibits at the higher energies clear discrepancies. The discrepancies between 
the theory based on pairwise forces only and data are largest in the region starting from the cross 
section minimum around Oc.m. ~ 130° up to 9c.m. ~ 180°. At energies up to about ~ 135 MeV 
these discrepancies can be removed when current three-nucleon forces, mostly of 27r-exchange 
character jsil, [3], are included in the nuclear Hamiltonian. At the higher energies, however, a 
significant part of the discrepancy remains and increases further with increasing energy. Especially 
complex picture exists for spin observables. Here adding 27r-exchange three-nucleon force into 
nonrelativistic calculations leads to eff^ects which depend on observable. They can be large or 
negligible, change their magnitude with energy and angle. Similarly to the elastic scattering cross 
section even after inclusion of three-nucleon force some of the discrepancies remain and increase 
with increasing energy. This indicates that additional three nucleon forces should be added to 
the 27r-exchange type forces. Natural candidates in the traditional meson-exchange picture are 
exchanges like n — p and p — p- This has to be expected since in xPT 39| in the order in which 



nonvanishing three-nucleon force's appear the first time there are three topologies of forces, the 
27r-exchange, a one-pion exchange between one nucleon and a two-nucleon contact interaction and 
a pure three nucleon contact interaction. They are of the same order and have to be kept together. 
Therefore it appears very worthwhile to pursue a strategy adding in the traditional meson exchange 
picture further three nucleon forces. Results presented here show that relativistic effects based on 
relativistic kinematics and boost effects of the nucleon-nucleon force play an important role in 
building up the magnitude of three-nucleon force effects. That gives hope, that taking the proper 
three-nucleon force into relativistic Faddeev calculations one will be able to improve the description 
of higher energy data for cross section and polarization observables. 
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Appendix A: Direct recalculation of the partial-wave projected three-nucleon force matrix 
elements from (p, q)- to (fc, q)- based basis 

We start from the matrix element of a three-nucleon force (p,q,a\v}^\l + P)\p' ,q' ju'^ in a 
partial wave basis used in norelativistic calculations with standard Jacobi momenta (p,q) [26] and 
would like to get the matrix element ^A:, g, a|V^^"'^^(l + P)\k' , q' , a''^ with p and p' replaced by the 
relative momenta of nucleons 2 and 3, k and k' , in their two-nucleon center of momentum system 



Using completeness of partial wave states one has 

'^k,q,a\v}^\l + P)\k',q,a'^ = J2 J dpp^ J dqfYl J dpP^ j dqq'^{k,q,a\p,q,a) 

a a' 

(p,q,d\v}^\l + P)\p,q,a) {p,q,a\k',q,a) . (Al) 



.2,q3j 



The partial wave state used in relativistic calculations |P, k, q, a) corresponding to the total three- 
nucleon center of momentum, P = 0, is given by 

\I>,k,q,a) = \P,k,q{l,s)j{X^)I{jI)JM;{t^)T)= E E E /rfq/dky;^,(k)iV(q; 

{^^,s\fi2,f^3,f^s)il,S,j\fii,fls,n) 

^i£2(^"''=(^-(-^)'^2(q2,q3))^if^3(^»c(5c(-g),/^3(q2,q3)) 

|q+ ip,W>|q2(k,-q) + ip,^'2)|q3(-k,-q) + ^P, /i'g)! (t, ^), T) (A2) 

where 

(q2,q3j = I p- I - 7^ r (A3j 

a[PNN k) UJff^^[PNN) Wfc Wfc 

is the Jacobian for the Lorentz transformation from (q2,q3) to (PjVAr,k) = (— q, k), = 



VmPT^, Mo = 2ujk = LOq^ + , and ojfj^ (Pnn) = y + Pnn ■ '^^^ momentum k2 (q2 , q3 ) = k 
and k3(q2,q3) = -k2(q2,q3). 

The nonrelativistic partial-wave state \P' ,p,q,a) with standard Jacobi momenta is given by 
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^,s\i>2,i>3,i>s){ls,j\i>i,i>s,i>){X, ^, i\i>x,i>i,i>i)U J, J\i>,Di,M) 
|q+ ip',z>i)|qr + ^P',^>2)|qr + ^P', ^^s)!^)^) (^4) 

where in the three-nucleon center of momentum system the nonrelativistic momenta q2"^ and ci^^ 
of the nucleons 2 and 3 are given by standard Jacobi momenta p and q as 

" TIT " ^ 

That leads to the scalar product (P, A;, Q;|P',p, g, a) 

(P,fc,g,a I P',p,^,a)= E E E /dq/^^ky4(k)iV(q2,q3) 

^ifM2 (^-c(i?c(-<?) , (q2 , q3))i^if;3 (i?«,c(i?c(-g), fes (q2 , qs)) 
Yxu.^ (q)(A, -,I\fix,fJ-i,Hi){j,I,J\n,Hi,M) 

E / / ^q^P,(P)^A,i>;,(q)(^'^'^l/^2,/^3>i^s)(M,i|^i,^s,^) 

2' '^l^'*" 

<5(q - q + ^(P - P'))'5(q2(k, -q) - ^d"' + ^(P " P')) 

(5(q3(-k, -q) - q3- + ^(P - P')) ^ (t^)r| (ti)r^ . (A6) 
That matrix element should be proportional to ^jj^^M independent from M. Thus 
(P,A;,g,a I P',p,g,d) = — ^E E E E f dci [ dk%^{k)N{ci2,cis) 

i^;£;(i?«,c(5e(-g), A;2(q2, q3))^if;3(i?«;c(i?c(-3), A;3(q2, qa)) 
5(q - q + ^(P - P0)5(q2(k, -q) - q2"'^ + ^(P - P')) 

<5(q3(-k, -q) - qr + ^(P - P'))((t^T)|(ti)f ) . (A7) 

The momenta of nucleons 2 and 3 in the three-nucleon center of momentum system are given 
through their two-nucleon center of momentum relative momentum k and the momentum of the 
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spectator nucleon 1, q, by 



and 



,3(-k.-<,) = -k-S__A_i_,. (AS, 
That allows to write the three (5-functions in the form 

S{q - q + 1(P - P'))«(q2(k, -q) - q^"' + ^(P - P'))i{q3(-k. -q) - + i(P - P')) = 
*(q - q + 1(P - P'))*(-P + k + ,^^p^^^^^^„^ q ^(P " P')) 

^(q-Wp-k- ^^^^^^^^^^^^ qWP-PO. (AlO) 

In the following we assume the three-nucleon center of momentum system (P = P' = 0) and 
drop the S(P — P') in all expressions. Performing the integration over dq one gets 

^(q2,q3)(2' 2^s\ii2, 11^, Hs){l, s, i\iiu lis, IJ) 
D'JiS^c{Bc{-q), fc2(q2, q3))i?if;3(i4;c(i?c(-9), fc3(q2, qs)) 

(q)(A, ^, /|/^A, /XI, /, Jl/x, /X,, M) 5^ y^,^ (q) 

^'^-■'- 2..(2!ir+Mo) ^><">l4>^'- 

That matrix element is a scalar which depends on the angles between the vectors q, k and p. These 
angles are fixed by the (5-function 5(p - k - ^-^|;^__q). Namely, while p = k + ^-^a_ci 
it follows that 

(k.q)2 



kp = kk + 
qp = qk + 



2ujk{2cuk + Mo) ' 

(q-q)(k-q) 

2cuk{2ujk + Mo) ' 



p p = p k+ (p q)^ . (ai2) 

Therefore one can take q pointing in z-dircction, what for given k, q, and p values, defines all 
angles between the appearing vectors. That allows to perform the integration over dq resulting in 

/I, I \ - g) ~ ^ 2a;fe(2tl+Mo)q|) 

{k,q,a\p,q,a) = —— / dx—- ^ '—^^ 2^ 2^ 2^ 



26 



-^(q2,q3)(^, ^, S\fi2, 1^3, l^s){l, S, jlfil, Us, 1^) 

D^J^ll{R^,{B,{-q), fc2(q2, q3))D'J^*^{R^,{B,{-q), h{q2, qa)) 

>A;,(q)(A,^,/|/^A,m,M/)(j,/,^l/^,M/,M) ^p,(p)^Az>,(q) 

i>si>ii/xi>i 

Q, i, J\D, Di, M){it, T\ {i, ^),f) (A13) 

with a; = q • k. We have chosen the coordinate system with q parallel to the z axis which leads to 
the components of q, k, p 

q = (0,0,g) , 



The isospin factor is 



k = {kVl-x^,0,kx) , 

p = (*yr^.O,te(l + ^;;^p|^)). (A14) 



it,h,T\ii,l-),T) = 5^i5j,f5MTM^ ■ ( A15) 



2 

Taking that all together gives 

-1 

(k)iV(q2,q3) 

^'Jil(^y^c{Bci-q), k2{q2, Ci3))D^J^;^{R^c{B,{-q), -fc3(q2, qa)) 

(2' 2'*l'"2''"^''"2 + /^3)(^si|M - m- H2- H3, 1^2 + 1^3, M - m) 

(A, ^, /|0, m, /, J|M - Ml, M)y;_^_^^_^, (^) 

(2 ' 2 ' ^1/^2' ^3, M2 + M3)(^~ j|M - //I - ^2 - M3> M2 + Ms, ^ " Ml) 

(A,^,/|0,Mi//i)(j,/,J|M-//i,//i,M) . (A16) 
The resulting expression for the matrix element {k,q,a\ v}^\l + P) \k,' q',a') is given by 



{k, q, a\vi^^ (1 + P)\k', q', a') = ^ 6,^6^ 



27rJ(2A + l)(2A + l) 



'^^ 2 J + 1 

a 

+1 

[dxYl E E^:M-Mi-M2-/.3(k)^^(q2,q3) 

^i-£*2 iRwc{Bc{-q), k2{q.2, ci3))Dy^*^{Ryjc{Bc{-q), A;3(q2, qs)) 

(2 ' 2 ' "'^1^2, Ms, M2 + M3)(^ s, - Mi - M2 - M3, M2 + M3, ^ - Mi) 
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(A, ^,I\0,ixi,ni){j,I,J\M - fii,ni,M)Yi^,j_^^_-^_-,^^{f>) 

2'^l'"2, A3,M2 + M3)(^,S,jW - Ml - A2 - A3,M2 +M3,M - ^l) 



~ 1 ~ . ~ „ 27rA/(2A' + l)(2A' + l) 
(A, 2> -'"lO, Ml, Mi)0', ^, - Ml, Ml, ^) 2^ ^t'i'^T'f' 2JTi 

a' 

+1 

(k')7V(q^,q'3) 

(i?^e(Sc(-9'), ^2 (q2, q3))^^£^ (i?^c(Sc(-g'), fc3(q2, Qs)) 
^, s'IM2, M3, M2 + M3)(^', s',j'\M -n'l- M2 - M3, M2 + M3, M - /x'l) 

(A', 2, -^'|0, Ml, Mi)(/, J\M - m'i, Ml, A^)^P,m-a*;-mJ,-a:,(p ) 

(2' 2' ^'IM2, Ms, M2 + M3)(^~', s',j'\M - m'i - M2 - Ms, M2 + Ms, M - m'i) 

(A', ^, 7'|0, m'i, M'i)(i', - Mi, m'i, M) 

^"-^ 2..(2!:r+Mo) -^ + 2..(2'l;iMo) -^^l'^'-^^) " ^^^^^ 

Appendix B: Transformation of a 3-dimensional three- nucleon force matrix element from 
(p, q) to (k, q) momenta 

We would like to express directly the matrix element (k, q | vj^'' |k', q') by that matrix element 
given in terms of single- nucleon momenta (qi,q2,q3| vj^'^ |qi,q2,q3). For momenta q2 and q3 of 
nucleons 2 and 3 in three-nucleon center of momentum system their relative momentum k(q2,qs) 
in the two-nucleon center of momentum subsystem of nucleons 2 and 3 is 



with Ei = \Jm? + q^"^. 

Using completeness of |qiq2q3) states one gets 



(k,q|14^^^|k',q') = j dqicZq2dqs5(qi + q2 + qs)(k, q|qi, q2, q3)(qi, q2, qslV'^ 
j (iqi(iq'2'iq'3|qiq2q3)(qi, q2, q'sik', q')^(qi + q'2 + qs) 
= j dqidq2dq3(5(qi + q2 + qs)(5(q - qi)(^(k - k(q2, qs)) 
j dqidq^dq^<5(qi + q^ + q^)<5(q' - qi),5(k' - k'(q'2,q'3)) 

= j dq2dq35(q + q2 + q3)<^(k - k(q2,q3)) 
j dq^dq^<5(q' + q^ + q'3)<5(k' - k'(q^, q'3)) 



a) 
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A^(q2,q3) W2,q3) 
= y" dk(q2, q3)(i(q2 + q3)^(q + q2 + q3)5(k - k(q2, qs)) 

J dk'(q'2, q'3)d(q'2 + q:j)<5(q' + q'2 + c^s^k' - k'{c^^, q'3)) 

iV(q2,q3)iV(q'2,q'3)(q,q2,q3ln^'^|q',q2,q3> • (B2) 

For given vectors k(q2,q3) = ko and q2 + qs = — qo the vectors q2 and qs are given by: q2 = 
— qo — q3 and qs = q!] with qg being the solution of the equation 

ko-k(-qo-q°,qO) = (B3) 

and similarily for the primed quantities. 
Thus one gets 

(kq|y«|k'q') = Ni-ci-4,ci)Ni-c{-4',4') 



q,-q-q3>q3 



q , -q - qs ' q3 



iV(-q-q^,q^)iV(-q'-q^',q^') 



1 

p = -2q-q3'q 



p' = -^q'-qr,q'^ ■ (B4) 



Starting from (jH]) and following the same steps as in ()B2p one gets for the partial wave projected 
matrix elements 

(A, , - , ^1 1/, /^/) (i, ,i,fii\j,fj-) ^ama (qi (^) 

1 ^ 



q2,q3) )] 

-^A3M3[^««=( Bc{-qi),h{ q2,q3) )] 

/ dk'dq'i ^ J2 XI (2''^2' o'/^3|s'>Ms)(^',/^hs',/^s>l/,/^j 

r-./ r-./ ../../ ^ ^ 



(A', /iA, -, A^'il/', M/)(/, ^/l-^, ^)^A'/.(,(qi)^'/.;(k') 

Se(-g'i),A:^(q'2,q3) )] 

Se(-g'i),A:3(q'2,q3) )] 
iV(q2,q3)iV(q2,q3)(q,q2,q3|Ki^'V',q2,q3) , (B5) 

where qi = gqi, k2(q2,q3) = A;k, k3(q2,qs) = -k2(q2,q3), and q2 together with qs result from 
()B3P and similarily for primed quantities. These partial wave matrix elements can be obtained 
using automatized partial wave expansion of Ref . [30|] . 
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FIG. 1: (color online) The matrix element of the TM99 3NF in relativistic- {{k,q,a\v}^^\k' ,q' ,a) - blue 
dashed line) and nonrelativistic-basis {{p, q, a\V^^^ \p' , q' , a) - red solid line) for the total angular momentum 
and parity of the 3N system J'" = and channel a = |(00)0(Oi)i(li)i >. The momenta p = k = 
0.132 fm-i and q = 0.132 fm-^. 
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FIG. 2: (color online) The same as in Fig. [T] but for the momenta p = fc = 5.25 fm ^ and q = 8.24 fm ^ . 
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2x10 




0246802468 10 



k [fm"^] k [fm"^] 

FIG. 3: (color online) The matrix element of the TM99 3NF in relativistic- ((k, qjl/j*-"'^'' |k, q) - blue dashed 
line) and nonrelativistic-basis ((p, q|1^4 |p, q) - red solid line). They are shown as a function of k at 
a number of q values assuming that all momenta are parallel to the unit vector ^ --i= , --i= , -1= ^ , all spin 
magnetic quantum numbers in the initial and final state are equal i, and isospins t = t' = Q. 
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FIG. 4: (color online) The elastic nd scattering angular distributions at the incoming neutron lab. energy 
E = 135 and 250 MeV. The solid (red) and dotted (blue) lines are results of the non-relativistic Faddeev 
calculations with the CD Bonn potential alone and combined with TM99 three-nucleon force, respectively. 
The relativistic predictions based on CD Bonn potential without Wigner spin rotations are shown by the 
dashed (blue) lines. The dashed-dotted (brown) lines show results of relativistic calculations with the TM99 
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three-nucleon force included. The pd data (x-es) at 135 MeV are from ref. f?| and at 250 MeV from ref. 
At 250 MeV also nd data of ref. t36i] are shown by circles. The inserts and figures in the right column display 
details of the cross sections in specific angular ranges. 
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FIG. 5: (color online) The vector (deuteron) Ay{d) and tensor analyzing powers Axx: ^yy^ and in 
elastic nd scattering at the incoming neutron lab. energy E = 70 MeV. For description of lines see Fig|31 
The pd data (open circles) are from [3]. 
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FIG. 6: (color online) The vector (deuteron) Ay{d) and tensor analyzing powers Axx: ^yy^ and in 
elastic nd scattering at the incoming neutron lab. energy E — 135 MeV. For description of lines see Fig|31 
The pd data (open circles) are from Q. 
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FIG. 7: (color online) The tensor analyzing powers -^xx -^yy -^zz 

in elastic nd scattering at the 

incoming neutron lab. energy E — 135 and E = 200 MeV. For description of lines see FigHl The pd data 
(open circles) are from [401] . 
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FIG. 8: (color online) The vector analyzing powers Ay{N) and Ay{d) in elastic nd scattering at the incoming 
neutron lab. energy E — 135 and E = 200 MeV. For description of lines see FigH) The pd data (open 
circles) are from 40 1. 
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FIG. 9: (color online) The spin correlation coefficients Cx.z and Cz^z and tensor analyzing power A^z in 
elastic nd scattering at the incoming neutron lab. energy E — 135 and 200 MeV. For description of lines 
see FiglH The pd data (open circles) are from 
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FIG. 10: (color online) The spin correlation coefficients Cx.x and Cz.x in elastic nd scattering at the incoming 
neutron lab. energy E = 135 and 200 MeV. For description of lines see FigH) The pd data (open circles) 
are from llOll. 
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FIG. 11: (color online) The spin correlation coefficients Czz,y and Cy^y in elastic nd scattering at the 
incoming neutron lab. energy E — 135 and 200 MeV. For description of lines see FigH) The pd data (open 
circles) are from [4^. 
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FIG. 12; (color online) The spin correlation coefficients Cxx,y — Cvv,v ^-nd Cxz,y in elastic nd scattering at 
the incoming neutron lab. energy E — 135 and 200 MeV. For description of lines see Fig|31 The pd data 
(open circles) are from |4o| . 
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FIG. 13: (color online) The spin correlation coefficients Cxy,x and Cyz,x in elastic nd scattering at the 
incoming neutron lab. energy E — 135 and 200 MeV. For description of lines see FigUl The pd data (open 
circles) are from [4^. 
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FIG. 14: (color online) The deuteron to neutron polarization transfer coefficients , i^Jf^, Kyy, and K^.^ 
in elastic nd scattering at the incoming neutron lab. energy E = 135 MeV. For description of lines see Fig|31 
The pd data (open circles) are from [35 1. 
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FIG. 15: (color online) The neutron analyzing power Ay{N) and neutron to neutron polarization transfer 
coefRcients , , , , and in elastic nd scattering at the incoming neutron lab. energy 
E — 250 MeV. For description of lines see FigU) The pd data (open circles) are from 137] . The nd data for 
Ay{N) (x-es) are from [sej] 
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FIG. 16: (color online) The five-fold cross section a / dnid^l2dE{°'^ for the breakup reaction d(n,np)n at 
El^^ = 200 MeV and fixed angles of out going nucleons 1 and 2 as indicated in the figures. For description 
of lines see FigH) The d(p,pn)p data (x-es) are from 41 {. 
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FIG. 17: (color online) The five-fold cross section (f" a / (Kli<Kl2dE^{^^ for the breakup reaction d(n,np)n at 
El^^ = 200 MeV and fixed angles of out going nucleons 1 and 2 as indicated in the figures. For description 
of lines see FiglH 
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FIG. 18: (color online) The polarization transfer coefficient Kjjy and deuteron analyzing powers Ay and Axx 
in the breakup reaction n(d,nn)p at incoming deuteron lab. energy E^j^^ — 270 MeV, shown as a function of 
the S-curve arc length. For description of lines see FigHl The 270 MeV dp data (circles) are from ref. 
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